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We reevaluate our published calculations of electromagnetic showers generated by GEANT 3.21 
and the radio frequency pulses they produce in ice. We are prompted by a recent report showing that 
GEANT 3.21-modeled showers are sensitive to internal settings in the electron tracking subroutine. 
We report the shower and pulse characteristics obtained with different settings of GEANT 3.21 
and with GEANT 4. The default setting of electron tracking in GEANT 3.21 we used in previous 
work speeds up the shower simulation at the cost of information near the end of the tracks. We 
find that settings tracking e~ and e + to lower energy yield a more accurate calculation, a more 
intense shower, and proportionately stronger radio pulses at low frequencies. At high frequencies 
the relation between shower tracking algorithm and pulse spectrum is more complex. We obtain 
radial distributions of shower particles and phase distributions of pulses from 100 GeV showers that 
are consistent with our published results. 
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In a recent paper |lj we reported the results of an ex- 
tensive study of the radio frequency pulse emitted by an 
electron-induced electromagnetic shower. Application to 
radio signal detection of ultrahigh energy electron neu- 
trino interaction in ice provided our motivation We 
used the standard particle physics detector-simulation 
package GEANT 3.21 (G3 afterwards) @ with default 
settings to generate the showers. Our total track-lengths 
of electrons and number of electrons at shower maximum 
were 68% of those in the paper by Zas, Halzen and Stanev 
(ZHS) Q|. Our pulse height at the Cherenkov angle at 1 
GHz was 72% of the ZHS result. Our study of sensitiv- 
ity of shower development to reasonable changes in input 
cross sections showed no effects large enough to explain 
the difference between G3 and ZHS shower profiles and 
track lengths. 

The source of the bulk of the discrepancy was recently 
presented by Alvarez-Muniz, Marquez, Vazquez and Zas 
5] . They noticed that the default setting of a parameter 
called "ABAN" in the electron tracking subroutine of G3, 
the setting used in our studies, gives shorter total, pro- 
jected and weighted track lengths and lower particle yield 
at shower maximum than GEANT 4 (G4 afterwards) or 
G3 with a choice of tracking rules different from the de- 
fault. Their results for track-length and shower profile 
from the latter simulations are 90% of those of the origi- 
nal ZHS work for these quantities and in closer agreement 
than G3 showers with the default setting Q, |{| . 

Track lengths are one indicator of the strength of elec- 
tric field pulses to be expected from the shower. However 
the electric field at a given angle and frequency depends 
on a combination of many factors. We make a compari- 
son among the shower and electric field pulse results from 
G3 (with several different combinations of electron track- 



ing rules) and G4. We report results for shower profiles, 
total track lengths, energy loss (dE/dx) and radial dis- 
tribution of the number of electrons minus the number 
of positrons around the shower maximum. We tabulate 
the G3, G4 and ZHS track-lengths and yields at shower 
maximum for ease of comparison of shower results. Re- 
garding the electric field in the Fraunhoffer regime, we 
report radio pulse strength at the Cherenkov angle and 
the phase behavior of the radio pulse both on and off 
the Cherenkov peak. We find that the particle yield and 
charge imbalance are larger with improved tracking. At 
low frequencies the pulses emitted at the Cherenkov an- 
gle are correspondingly stronger as well. 

Electron tracking in GEANT: The subroutine respon- 
sible for tracking electrons in G3 is called GTELEC. In 
its default setting, ABAN=1, if the electron range is too 
short to reach the boundary of a detector volume, and the 
distance to its next bremsstrahlung emission is greater 
than its range, the electron's energy is deposited in the 
volume and it is stopped. This is done regardless of the 
value of the electron kinetic energy threshold parame- 
ter CUTELE set by a user in the input card file. This 
was presumably made the default to increase the cal- 
culation Al speed of high energy particle physics detec- 
tor simulations in non-sensitive volumes, where detailed 
tracking was not needed 6J. With the default setting, 
however, track lengths are reduced significantly, with a 
consequent reduction in the calculated radio Cherenkov 
emission. This is what we observed in our previous study 
Q , where we compared results from G3 with default set- 
ting to the results of 0. 

In practice, it seems for applications such as ours the 
preferred setting is ABAN=0, which tracks all elec- 
trons down to the threshold set by the user. The setting 
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ABAN=2 allows the user to track electrons precisely (sim- 
ilar to ABAN=0 setting) in sensitive volumes or not so 
precisely in the non-sensitive volumes 0. "Sensitive" vs. 
"non-sensitive" settings are distinguished by the value of 
a step parameter called "STMIN" in the tracking medium 
parameter list, where "sensitive" produces many more 
steps per track than "non-sensitive" . G4 has range cuts 
that correspond to energy thresholds for the production 
of particles. A particle, once produced, is tracked down 
to rest regardless of its production threshold. 

Results: We expand on the study reported in Ref. 
where showers produced with G3 ABAN=1, G3 ABAN=2 
with "sensitive" volume, G4 and ZHS are compared. 
We run G3 for four different tracking choices, namely 
ABAN=0, ABAN=1, ABAN=2 sensitive volume (S) and 
ABAN=2 non-sensitive volume (NS) settings, to create 
showers in ice. We run G4 for comparison as well. We 
then repeat our study of the key features of the show- 
ers and the consequent key features of the radio pulses: 
pulse angular distributions, phase distributions and fre- 
quency spectra. 

1 ) Shower Characteristics: First we investigate the im- 
portant issue of total track lengths in Table [I] We note 
that the track lengths for ABAN=2 (NS) setting are sim- 
ilar to the default ABAN=1 setting; the ABAN=2 (S) set- 
ting used in Ref. [f| gives results similar to the preferred 
ABAN=0 setting and both are similar to the results from 
G4. The agreement between track lengths produced by 
these latter three GEANT versions and those produced 
by the ZHS code improves as one compares total, total 
projected e + p and total projected e — p track lengths. 
Comparing the projected e+p and e—p track lengths, we 
see that the loss of track length for ABAN=1 is more sig- 
nificant for electrons than for positrons. We expect this 
to be the case, since a large part of the electron excess 
comes from low energy electrons swept into the shower 
by Compton scattering and S-iay production. 

Despite the improved agreement on total track lengths, 
there are still significant differences between the G3 
ABAN=0, G3 ABAN=2(S), and G4 simulations. For each 
code individually, however, the variance in total track 
length 36 — 40 cm is reasonable. For a given shower en- 
ergy Eq-, the total track length L to t is determined pri- 
marily by the rate of ionization energy losses, where the 
energy transfer is less than a low threshold E t ~ 100 keV. 
Then L to t ~ Eo/(dE/dx)i on , as expected under the con- 
ditions that a) soft ionization losses are essentially in- 
dependent of particle energy and b) hard interactions do 
not "remove" energy from the shower. Hard 5-rays create 
their own track length and hard 7's return their energy 
to the shower in pair production and Compton scatter- 
ing processes. Some energy is lost from the shower as soft 
particles fall below E t . The fluctuations in track length 
are thus determined by the fluctuations in the number 
of soft particles multiplied by E t /(dE/dx) ion . Estimat- 
ing the number of soft particles by N t = Eo/E t , the 
fluctuations in total track length are of order 8L to t ~ 
E t /(dE/dx), on y/W t = ^/E E t /(dE/dx) lon = 50 cm, for 



E = 100 GeV, E t = 100 keV, and {dE/dx) ion = 
2 MeV/cm. Although this is a rough estimate, the ob- 
served values of a for L to t are consistent with expecta- 
tions for the newer simulations. 

This estimate leads to several observations, a) The 
deviations in L tot between G3 ABAN=0, G3 ABAN=2(S), 
and G4 are small but significant; this is under study, 
b) Due to the default decision in the tracking code to 
abandon some particles, the older G3 simulations have 
an effective value of Et ~ 5 MeV, which implies 8L ~ 
3.5 m, generally consistent with the Table entries, c) 
For e — p track length, the fluctuations are dominated by 
the competition between pair production and Compton 
scattering, for which one may estimate E t ~ 20 MeV and 
5L ~ 7 m. 

Particle yields at the maxima of longitudinal shower 
profiles are also listed in Table [I] The results from the 
ZHS code, G3 with ABAN=0, ABAN=2 with sensitive 
volume, and G4 all agree reasonably well. Agreement is 
especially good for the excess charge, where results agree 
within errors. Again, the low yields from G3 ABAN=1 
and ABAN=2 (NS) are consistent with the interpretation 
that lost track length is caused by a deficiency of low 
energy electrons being scattered into the shower. 

Next we consider how individual particle tracks con- 
tribute to the radio pulse. For low frequencies it can 
be shown that the contributions are in phase and the 
pulse is proportional to the weighted, projected track 
length, labeled (e-p) in Table 1. As the frequency in- 
creases, this simple coherence breaks down as described 
by two parameters, which we called 1] the Cherenkov 
Phase (CP) and the Translational Phase (TP). The CP 
is related to the phase coherence within a single track seg- 
ment, whereas the TP describes the relative phase among 
track segments due to their different positions within the 
shower, relative to the observer. 

In general, we find that the different simulations pro- 
duce phase distributions similar to those shown in Fig. 
17 of Ref. lj. Specifically, the radial charge distribution 
at shower maximum, shown to peak at p ~ 0.5 cm in 
Fig. 22 of Ref. 0, is reproduced apart from an overall 
increase in normalization 9]. Moreover, the radial charge 
distribution early in the shower evolution has a core dom- 
inated by high energy particles, and the normalization of 
this core does not depend on the details of the electron 
tracking at low energy. 

2) Radio frequency pulse: Next we examine how the 
changes in shower simulation affect the modeled radio 
Cherenkov pulses. At low frequencies, all particles con- 
tribute in phase, so the increase in the excess track length 
found in G3 ABAN=0 and G4 should be reflected by an 
increase in the amplitudes of the low frequency part of 
the spectrum as compared to showers produced with G3 
ABAN=1. In Fig. \I\we see that this is indeed the case, 
where the spectrum R x \E\ is shown for an observer on 
the Cherenkov cone. 

In Fig|3we extend the frequency scale out to 10 GHz, 
which suffices to show clearly a transition from low to 
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TABLE I: Track length and particle yield results from an average of 100 GeV electron induced shower using different Monte 
Carlo shower codes. The error bars correspond to error in the mean, s/y/N, where s is the standard deviation reported within 
parentheses and TV is the number of showers used in each case. 



Monte Carlo 
shower code 



No. of 

showers 



G3 (ABAN=0) a 20 

G3 (ABAN^l) 6 20 

G3 (ABAN=2, NS) 20 

G3 (ABAN=2, S) c 20 

G4 100 

ZHS 20 



Absolute 
(e + p) [m] 



Total track lengths 
Projected 
(e + p) [m] 



Projected 
fe- 



rn 



Particle yield at shower max 
Total Excess Excess 

(e + p) (e — p) fraction 



542.74 
389.51 
389.40 
561.05 
572.58 
642.17 



0.08 (0.36) 
0.48 (2.16) 
0.71 (3.18) 
0.09 (0.40) 
0.04 (0.40) 
0.06 (0.25) 



455.3 ± 

365.7 ± 
367.0 ± 

459.8 ± 
466.3 ± 
516.6 ± 



0.2 (0.9) 
0.5 (2.2) 
0.7 (3.1) 
0.3 (1.3) 
0.2 (1.6) 
0.2 (0.9) 



125.0 ± 2.0 (9.0) 

76.3 ± 1.5 (6.7) 

78.4 ± 2.8 (12.5) 
125.8 ± 1.9 (8.5) 
135.0 ± 0.8 (7.9) 
135.2 ± 1.5 (6.8) 



148 ± 5 (22) 42 ± 3 (13) 28% 

111 ± 7 (31) 20 ± 2 (9) 18% 

113 ± 10 (45) 25 ± 5 (22) 22% 

148 ± 8 (36) 41 ± 3 (13) 28% 

153 ± 3 (27) 45 ± 1 (13) 30% 

164 ± 6 (28) 44 ± 2 (11) 27% 



"Preferred setting (new). 
'Default setting (old). 

c Setting used in Ref. H (similar to ABAN=0). 



high frequency for the same sets of simulations. In this 
range, the spectra for the G3 ABAN=0 and G4 showers 
flatten out, but the G3 ABAN=1 spectrum grows slowly, 
so that the difference in their normalizations, relative to 
G4, is less than 20% for / = 10 GHz. We infer that in 
this region the spectra are not simply described by the 
phase coherent sum of amplitudes from a large number of 
low energy electrons. We note two effects that could be 
important in this region. First, a relatively small number 
of high energy particles near the shower axis could con- 
tribute strongly to the spectrum in this region. Since this 
part of the shower evolution is not sensitive to low en- 
ergy electron tracking, we may see the G3 and G4 results 
converging. This trend continues at higher frequency. 
Second, a large fraction of the track-segment amplitudes 
are no longer coherent. For that fraction, the magnitude 
of E increases as a random walk as the track length in- 
creases. In this case, one would expect the amplitude 
to increase as the square root of the number of track 
segments, rather than linearly, which reduces the ratio 
of the pulse from the G4 generated shower to that from 
ABAN=1 as the frequency increases. This is seen clearly 
in Figs. 1 and 2. Given that total and projected track 
lengths are proportional to energy 0,3, the effects just 
outlined should be relevant when comparing pulses from 
showers of different energies using the same simulation. 

Figure [3] shows a comparison between pulse spectra 
from 1 TeV showers and 100 GeV showers, both gener- 
ated by G4. The spectra are extended up to 50 GHz. For 
frequencies above 15 GHz, the results are fully consistent 
with Fig. 25 of Ref. lj. At low frequencies, below a few 
GHz, the results are consistent as well, after account- 
ing for the change in normalization between G4 and G3 
ABAN=1. At 50 GHz, the TeV average, scaled down by a 
factor of 10, is smaller by roughly a factor 3 ~ \/l0 than 
the 100 GeV average, possibly hinting that the "random 
walk" increase of E with energy, and therefore total track 
length, is setting in. At the same time, it seems clear that 
a combination of effects must be in play for intermediate 
frequencies. 




Frequency (MHz) 

FIG. 1: Direct Monte Carlo calculation of the low frequency 
end of the spectrum of the average electric field amplitude 
for twenty 100 GeV showers. The showers are run using G4 
and G3 with both the preferred ABAN=0 and the default 
ABAN=1 settings. The same code was used to calculate elec- 
tric fields from all the GEANT showers. 



Since the flat spectra observed above several GHz con- 
trast with the sharp falloff expected above the 1-2 GHz 
scale 0, , we referred to this phenomenon as an "ex- 
tended coherence regime" in [1]. This is different from 
the use of coherence to refer to the rapid rise of the elec- 
tric field strength as frequency increases in the region 
where frequency is low compared to the scale set by the 
peak at 0.5 cm in the shower radial distribution. 

Summarizing, we see that the improved electron track- 
ing performed by G4 or G3 ABAN=0 results in propor- 
tionately increased radio Cherenkov spectra at frequen- 
cies below a few GHz. At high frequencies the situation 
is more complicated. We are conducting further studies 
at a variety of energies and including showers initiated 
by 7-rays to identify the determining factors. 

Summary and further discussion of results: 

Motivated by the observation in [f| that the default 
setting for electron tracking in G3 stops many electrons 
prematurely, we have reviewed the calculations reported 
inQ. In pj we used G3 with default setting and found 
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FIG. 2: Direct Monte Carlo calculation of the frequency spec- 
trum of the average electric field amplitude for twenty 100 
GeV showers. The showers are run by G4 and G3 with both 
the preferred ABAN=0 and the default ABAN=1 settings. 
The same code was used to calculate electric fields for all the 
GEANT showers. 
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FIG. 3: Direct Monte Carlo calculation of the frequency spec- 
trum of the average electric field magnitude for one hundred 
100 GeV showers and ten 1 TeV showers run by G4. The 1 
TeV spectrum is divided by a scaling factor of 10. 



a persistently smaller shower and a weaker electromag- 
netic pulse at low frequencies than those reported in [4J. 
We find that this discrepancy largely disappears when 
electrons are tracked more completely in G3 or when us- 



ing G4 generated showers, similar to results reported in 
13 ■ Figure 1 shows that the strength of the pulses at low 
frequency clearly reflect the track lengths and particle 
yield at shower maximum. The high frequency behavior 
is more complex and is under further study. 

The particles in the electromagnetic showers in a uni- 
form medium, ice in the case of interest here, should be 
tracked as long as they still contribute to the Cherenkov 
signal. This is the case, for the ZHS, G4 and G3 sim- 
ulations with the tracking option ABAN=0 or ABAN=2 
with sensitive volume; the results from all these simula- 
tions agree at the 10% level. There is physics, namely 
significant emission, that is lost if the tracks are stopped 
prematurely and their energy dumped into the medium. 
Though choosing the default (ABAN=l) setting affects 
only the low energy tail of the shower, below the few 
MeV range, the population of particles (predominately 
electrons) is so high at the low energy end that the field 
contribution cannot be ignored. 

For applications of the shower and pulse codes to neu- 
trino detection experiments, the difference is obviously 
important. For example, below 1 GHz where currently 
deployed antennas are sensitive, the expected radio emis- 
sion pulse is 1.3 to 1.4 times higher using the G3 or G4 
code with particles tracked to threshold compared to the 
truncated showers produced by G3 with the ABAN=1 
setting. The more correct treatments that follow the par- 
ticles to lower energy produce significantly larger fields 
on the Cherenkov cone. This increases the effective vol- 
ume estimates and leads to higher sensitivity estimates 
for cosmic neutrino detection systems like the Radio Ice 
Cherenkov Experiment (RICE) ^l| and the recently ap- 
proved balloon-borne experiment ANITA , both using 
the Antarctic polar ice cap as a neutrino target. 
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